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1.  Background 


A  complex  mechanical  system  analysis  can  incur  several  hundred,  thousand,  or  even  million 
degrees  of  freedom  (DOF).  As  the  number  of  DOF  in  the  system  increases,  so  does  the 
computing  time,  in  some  cases  exponentially.  Consequently,  computational  engineers  are 
always  seeking  to  strike  a  balance  between  model  accuracy  and  computational  efficiency. 

In  1960,  Hurty  (1)  introduced  the  concept  of  modal  synthesis,  also  called  substructure  or 
component  mode  synthesis.  It  has  enjoyed  popularity  over  four  decades  due  to  its 
straightforward  methodology.  The  premise  of  modal  synthesis  is  that  the  first  few  modes  of 
vibration  of  a  complex  structure  can  be  predicted  by  using  the  lower  modes  of  vibration  of  its 
substructures  (2).  An  extension  of  modal  analysis,  modal  synthesis  provides  for  a  reduction  in 
the  overall  number  of  DOF.  As  such,  the  computational  complexity  decreases,  but  the  reduced 
model  still  maintains  the  desired  numerical  precision. 

The  theory  behind  modal  analysis  is  fairly  simple.  A  coordinate  transformation  of  the  equations 
of  motion  yields  an  equivalent  uncoupled  system  of  equations,  each  of  which  is  linearly 
independent  of  the  others.  The  transformation  yields  the  equations  in  terms  of  principal 
coordinates  rather  than  physical  coordinates.  The  physical  motion  of  the  system  then  is  simply  a 
linear  combination  of  the  motion  of  the  principal  coordinates  such  as  in  equation  1 . 

=  [(p(2;)][u(t)] 

{q{z,t)}  =  +  (?,{z)v,{t)  +  ...  +  (p^{z)v^{t) 

where  v,(t)  is  the  principal  coordinate  response  of  the  diagonalized  equations  of  motion,  and 
{(p;(z)}  is  the  corresponding  eigenvector  or  mode  shape. 

One  benefit  of  modal  analysis  is  that  it  shows  the  contribution  of  the  individual  modes  of 
vibration  to  the  motion  of  the  entire  structure  and  thus  explicitly  captures  the  dynamics  of  the 
system.  The  lower  modes  extracted  by  modal  analysis  generally  provide  the  most  contribution  to 
the  motion  of  the  system.  Therefore,  it  provides  a  basis  for  DOF  reduction  and  model 
simplification.  Using  modal  analysis,  the  following  applications  are  available:  among  others, 
structural  modification,  sensitivity  analysis,  model  reduction,  correlation  of  model  and  finite 
element  results,  response  prediction,  substructure  coupling,  and  structural  damage  detection  (5). 

In  dynamic  analysis,  traditional  modal  analysis  techniques  have  achieved  varying  degrees  of 
success  when  a  system  undergoes  damping,  due  to  the  presence  of  the  damping  matrix.  In 
structural  applications,  the  common  approach  is  to  assume  proportional  damping,  such  that 
traditional  eigenvalue  approaches  still  apply.  In  rotor  dynamic  applications,  however,  the 
damping  is  rarely,  if  ever,  proportional. 
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Furthermore,  the  gyroscopics  of  a  spinning  rotor  produee  frequeney  variations  which  are  absent 
in  non-spinning  systems.  One  of  the  first  rotor  dynamic  analyses  to  use  modal  synthesis  resulted 
in  the  computer  program  Analysis  of  Rotor  Dynamic  Systems  (ARDS  )  that  could  return  the 
frequency  modes  within  1-10%  error  {4).  However,  gyroscopic  effects  were  not  explicitly 
defined  within  the  program. 

While  the  gyroscopic  forces  are  themselves  conservative  and  take  no  energy  from  the  system, 
they  can  become  significant  at  higher  rotational  speeds  and  must  be  addressed. 

The  presence  of  gearing  also  complicates  the  modal  analysis  process.  Separately,  the  scientific 
studies  of  gearing  and  rotor  dynamics  first  appeared  approximately  150  years  at  the  height  of  the 
Industrial  Revolution.  Studied  together,  geared  rotor  dynamic  models  did  not  appear  until  the 
1960s  (5).  Gears  couple  the  motion  of  rotating  shafts  together  and  result  in  coupled  lateral- 
torsional  or  lateral-torsional-axial  vibrations  not  found  in  non-geared  rotating  shafts. 

Additionally,  the  gear  mesh  induces  a  periodic  forcing  function  known  as  the  gear-mesh,  or 
tooth-pass,  frequency.  The  gear-mesh  frequency  has  a  unique  impact  and  is  further  discussed 
later. 

Returning  to  the  subject  of  general  damping  in  dynamic  models,  previous  attempts  have  made 
approximations  by  different  methods  (6),  the  most  common  being  the  proportional  damping 
assumption.  Once  again,  this  is  not  particularly  useful  in  rotor  dynamics.  Another  approach 
ignores  damping  in  the  system  and  focuses  only  on  the  mass  and  stiffness  contributions,  resulting 
in  the  familiar  eigenvalue  problem:  (K  -  M  =  0).  A  third  method  has  been  to  transform  the 
damping  matrix  via  the  traditional  modal  matrix  and  simply  ignore  the  off-diagonal  terms.  Other 
methods  (7,  8)  perform  an  iterative  procedure  in  transforming  systems  consisting  of  symmetric 
damping  matrices  only.  However,  a  modal  transformation  does  exist  (9,  10)  to  uncouple  the 
equations  of  motion  for  a  system  undergoing  general  damping  that  can  include  gyroscopic 
effects.  This  transformation  procedure  provides  the  theoretical  basis  for  the  proposed 
methodology. 

The  overarching  objective  of  this  report  is  to  present  a  method  for  conducting  modal  synthesis  on 
a  geared  rotor  dynamic  system  under  the  influences  of  non-proportional  damping  and  gyroscopic 
effects.  In  addition  to  model  simplification,  this  technique  also  exhibits  further  potential  for  use 
in  optimization  and  system  identification  schemes. 
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2.  Coordinate  Transformation 


2,1  Theory 

For  a  rotor  dynamic  system  undergoing  general  damping  and  gyroscopie  effeets,  the  general 
equations  of  motion  are 

lMl{q}  +  [D]{q}  +  lK]{q}={F(J)} 

|D|  =  |C]  +  $2[G| 

where  [D]  is  the  “damping”  matrix,  eonsisting  of  the  general  damping  matrix,  [C],  and  the  skew- 
symmetrie  gyroseopie  matrix,  [G].  The  gyroseopie  effeets  in  the  system  are  a  funetion  of  the 
shaft  rotational  speed,  (f2),  and  the  gyroseopie  matrix.  This  representation  is  often  defined  as  the 
“real”  representation,  whieh  will  double  the  size  of  the  system.  Furthermore,  distinction  between 
forward  and  baekward  whirl  disappears  using  this  representation  {11). 


The  presenee  of  the  damping/gyroseopie  matrix,  [D],  poses  a  problem  for  traditional  analysis, 
sinee  the  eigen-solution  problem  uses  only  two  matriees.  Converting  the  problem  into  the 
familiar  first-order,  state-spaee  equation,  the  general  eigenvalue  problem  beeomes 


T  OUql  [0  lUql  [0^ 
0  -K  |q|  ^  M  D  q  “  0 


(3) 


or 


where 


-  Aj{q'}  +  [A2]{q'}  =  {0} 


0  I 
M  D 


Assuming  the  displaeement  and  veloeity  veetors  in  equation  4  to  be  of  the  form; 

{q'}  =  {Qje" 

{q'}  =  A{Q,}e" 


the  generalized  eigenvalue  problem  for  the  state-spaee  formulation  beeomes 

([Aj-A[A,]){(t,J  =  {0} 


(4) 

(5) 

(6) 

(V) 
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where  (A)  is  the  eomplex  eigenvalue,  and  is  the  eorresponding  eigenveetor.  The  state- 
space  formulation  results  in  (2N)  eigenvalues  occurring  in  complex  conjugate  pairs,  where  (N)  is 
the  number  of  DOF  of  the  system. 

Returning  to  equation  2,  the  system  matrices,  [M],  [K],  and  [D],  are  generally  not  diagonal  but 
have  cross-coupling  terms  throughout.  The  objective  is  to  determine  the  mass,  stiffness,  and 
damping  matrices,  respectively  as  [M/j],  [Kz)],  [Dzi],  that  are  diagonal  matrices  of  the 
transformed  equations  of  motion. 


The  next  step  is  to  calculate  the  eigenvalues  and  eigenvectors  of  the  state-space  formulation. 
The  results  will  form  two  modal  matrices,  m  and  ['F/j,  whose  columns  contain  the  right  and 
left  eigenvectors,  respectively. 

Since  the  eigenvalues  occur  in  conjugate  pairs:  A,  =  a± jP,  the  eigenvalues  and  corresponding 
eigenvectors  can  arrange  so  that  the  conjugate  pairs  are  separate.  From  this,  two  diagonal 
matrices  of  dimension  (N  x  N)  can  be  defined  that  contain  the  conjugate  eigenvalues  on  the 
diagonal: 


A 


1 


ctj+j/Jj  0  0 

0  ■•.  0 

0  0 


A 


2 


0  0 
0  ■•.  0 

0  0  %-j/3^ 


(8) 


(9) 


If  any  of  the  eigenvalues  contain  purely  real  roots,  the  transformation  scheme  changes  somewhat 
(P).  However,  for  brevity,  all  roots  are  assumed  to  be  complex. 


To  ensure  the  desired  results,  the  right  eigenvector  must  be  scaled  or  normalized.  One 
appropriate  scaling  results  in 


%  AJIT 


(10) 


such  that 


T 

A, 

0 

'P, 

A, 

= 

= 

1 

r 

0 

A. 

T 

I  0 

% 

A  ' 

= 

= 

I 

1 

2 

r 

0  I 

(11) 
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Both  the  left  and  right  eigenvectors  are  complex,  and  any  further  transformations  will  continue  to 
produce  complex  numbers.  The  elimination  matrix,  [S],  is  therefore  defined  to  eliminate  all 
complex  notation  from  the  transformation. 


—  I  ^ 

42  42 


(12) 


Then 


T 

T 

L 

L 

3 

a/ 

3 

= 

X 

y 

1 

1 

r 

L 

L 

y 

z 

T 

T 

I  0 

M 

A  ' 

M 

1 

2 

r 

0  -I 

(13) 

(14) 


The  matrices  L 


and  L  are  real  and  diagonal.  They  are  equal  to  the  following: 


=  rea/(Aj 
=  imag{K^) 

L  =  —real(A,)  —  —A 

2  ''1^  X 


A  diagonal  matrix,  [F],  is  defined. 

T{k,k) 


2"'''  2L^(A:,A:) 


Then  the  following  relationship  holds. 


cosh  r  sinh  F 

L  L 

X  y 

cosh  F  sinh  F 

0 

CO 

n 

sinh  F  cosh  F 

L  L 

y  z 

sinh  F  cosh  F 

CO 

n 

2C« 

% 

(15) 


(16) 


(17) 


where  the  matrices  [co„]  and  [2(!^co„]  are  diagonal,  containing  the  natural  and  damped  natural 
frequencies  of  the  system,  respectively. 


The  left  and  right  transformation  matrices  can  now  be  determined.  These  are  the  matrices  that 
diagonalize  the  state-space  matrices  [Ai]  and  [A2]. 


coshF 

sinhF 

CO 

L  ^ 

0 

sinhF 

coshF 

0 

I 

(18) 
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cosh  r  sinh  F 

CO 

0 

M 

n 

r 

sinh  r  cosh  F 

0 

I 

(19) 


Transforming,  the  diagonalized  results  become 


t/ 

T 

A,' 

T  ' 

B  ' 

0  K, 

L 

1 

R 

1 

K,  D, 

T 

0 

t/ 

A  ' 

T  ' 

= 

B 

= 

D 

2 

R 

2 

0 

-M, 

(20) 


(21) 


This  transformation  methodology  will  result  in  {2N)  de-coupled  linearly  independent  equations. 
The  usual  relationships  achieved  through  modal  analysis  apply  here  as  well,  namely 


M, 


I 


K.  =  K] 


(22) 


If  carried  through  from  start  to  finish,  the  process  produces  results  consistent  with  the  following 
simplification  {12). 


\  =  a  +  ](3 

4  j/3 


^i-~  i^i  +  ^2 )  -  ~‘^^i 

^  _  (^2  +  ^l)  ~  (^2  ~  ^l) _ 2 


m.  =  1 


(23) 


where  i=\  ...  N 

From  modal  analysis,  the  transformed  coordinates,  {v},  are  principal  coordinates.  The 
transformation  to  the  physical  coordinate  system  is 


q 

T  ' 

V 

q 

R 

u 

(24) 


In  the  equation  above,  the  term,  {u}  is  not  necessarily  the  time  derivative  of  the  general 
coordinate,  {v},  except  when  the  system  is  in  free  vibration  (9).  Also,  the  vector  order  of  the 
principal  coordinate  and  its  “derivative”  are  reversed  compared  to  the  physical  coordinates. 
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2.2  Transformation  Simplification 

The  full  procedure  outlined  in  the  previous  section  is  rigorous  and  lengthy.  However,  a 
simplified  transformation  is  available  using  equation  23  and  similarity  transformation 
relationships.  From  the  definition  of  a  similarity  transformation, 

(25) 

'*'J"[A2l['p:l=[il=[®zr[B.l[®:l  (26) 

where  m  and  [O/]  are  the  modal  matrices  of  the  reduced  system.  Using  linear  algebra,  the 
transformation  matrices  [T^]  and  [Tj;]  can  be  written  as  (12) 

T,]  =  [^.]Kr  (27) 

Tl  =  [T'l[o'r  (28) 

Using  these  relationships,  the  simplified  transformation  procedure  is  as  follows: 

1 .  Given  matrices  [Ai]  and  [A2],  solve  the  eigenvalue  problem.  Arrange  the  eigenvalues  in 
matrices  [Ai]  and  [A2].  Assemble  the  left,  ['?/],  and  right  normalized,  'P^  ,  eigenvector 

matrices,  corresponding  to  [Ai]  and  [A2]. 

2.  Using  equation  23,  calculate  the  diagonal  values  of  [Mz)],  [K/)],  [Dz)].  Substitute  into 
equations  20  and  21  for  [Bi]  and  [B2],  respectively. 

3 .  Determine  the  left,  m,  and  right,  O  I ,  normalized  eigenvector  matrices  of  the  modal 
system,  [Bi]  and  [B2]. 

4.  Substitute  into  equations  27  and  28  to  solve  for  [Tz]  and  [T^],  respectively. 

In  order  to  ensure  more  precise  results,  it  is  possible  to  further  refine  the  transformation  by 
calculating  the  left  eigenvector  matrix,  [T'/],  using  either  equation  20  or  21,  rather  than  using  the 
originally  determined  left  eigenvector  matrix.  The  results  are  then  substituted  into  equations  27 
and  28. 
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3.  Modal  Synthesis 


3.1  Theory 


In  modal  synthesis,  the  analysis  begins  with  the  substruetures.  The  equations  of  motion  for  eaeh 
substrueture  are 


,1K}+[d,|K}+[k,|K}  =  {f(o.} 


where  the  subseript,  (z),  represents  the  z*  substrueture.  The  substrueture  eigenvalue  problem 
solution  follows  the  state-spaee  methodology  introdueed  in  equation  3. 

I  0  I  [q|  [O  l|  |q|  [o’ 

~  0  -K  q  ^  M  D  q  “  F(t)  ' 


Onee  the  substrueture  analysis  is  eomplete,  the  global,  uneoupled  system  is  eonstrueted. 
However,  unlike  some  other  modal  synthesis  works  {13),  the  entire  state-spaee  formulation  of 
equation  30  beeomes  part  of  the  global  model.  The  dimensions  of  the  global  state-spaee 
formulation  are  [(2M  +  . . .  -H  INm)  x  {2N\  +  . . .  +  2Nm)\ 


T  0 
0  -K 


0  I 
M  D 


T  0 
0  -K 


0  I 
M  D 


or  in  simpler  notation 


“  I  W).  [[Asl,  0  I  {qf 


{F(t)}, 


1  »  lA.Llllq'll  1  0  1^ 

where  (m)  is  the  number  of  substruetures  in  the  system. 


{q’L 


For  simplieity,  equation  32  ean  be  rewritten  as 


A:]{q}  +  [A;]{q}  =  {F(t)‘) 
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where  the  asterisk  symbolizes  the  global  system  of  individual  substruetures.  Note  that  the 
substrueture  eoupling  has  not  been  introdueed  yet. 


Returning  to  the  fundamental  eoncept  of  modal  analysis,  the  system  equations  of  motion  are 
transformed  into  principal  coordinates.  The  global  transformation  matrices  are  composed  of  the 
substructure  transformation  matrices  of  equations  27  and  28,  respectively. 


f  *  1 
jv  > 


+ 


r  ♦  *  1 
y 


(34) 


where 


(35) 


Equation  34  simplifies  to 


B 


{fW*} 


(36) 


or  in  long  form, 
0  K„ 


0  K, 


K„  0 
0  -M 


K,  0 

0  -M 


=  T. 


0 

F(t)j 

F(t) 


(37) 


Note  that  the  global  equation  containing  the  (m)  systems  is  still  uncoupled. 

The  coupling  occurs  as  a  result  of  the  gear  mesh.  The  mesh  stiffness  is  the  coupling  agent. 
However,  this  treatment  is  general  enough  to  apply  to  other  types  of  couplings  with  only  minor 
changes  to  the  methodology. 

The  mesh  matrix,  \K.mesh\,  breaks  into  quadrants  as  shown  below.  A  methodology  for 
determining  the  components  of  the  mesh  matrix  for  spur  or  helical  gears  can  be  found  in  {14). 

|K„„„ 


[k 

lii 

[k 

[  mesh  1 

[  mesh 

[k  ] 

ji  1 

K  ] 

[  mesh  \ 

mesh  J 

(38) 


For  a  two-shaft  system,  the  coupling  becomes 


9 


I  0 
0  -K 

0 


0 

I  0 
0  -K 


-  K 


mesh 


-  K 


mesh 


-  K 


mesh 


hJ 


-  K 


mesh 


(39) 


or 


+ 


mesh 


(40) 


Using  the  same  transformation  on  the  seeond  term  in  equation  39  results  in  the  transformed 
eoupling  matrix,  \KDm\- 


B, 


mesh 


]0 

0 

0 

0 

0 

0 

0 

0  ]' 

0 

[^mesA 

0 

_ 1 

L  TTiesh  j 

0 

0 

[0 

0 

0 

0 

0 

0 

0 

0 

0 

_ 1 

L  mesh  J  j 

0 

— [k  1 

L  mesh  J  J.  J 

0 

0 

(41) 


Note  that  sinee  this  is  a  eoupling  matrix,  it  is  not  neeessarily  diagonal.  Thus  the  eoupled  system 


modal  matrix. 


B 


IS 


B 


sys 


B 


+ 


B 


mesh 


(42) 


In  this  formulation,  it  is  assumed  that  no  damping  oecurs  in  the  gear  mesh.  Consequently,  the 
gear  mesh  does  not  affeet  the  modal  system  matrix,  B*  from  equation  36.  It  requires  no 

modification: 


B„ 


B 


(43) 


Equation  36  then  becomes 


(44) 


The  general  eigenvalue  problem  of  equation  44  yields  the  same  results  as  the  original  problem  of 
equation  3. 
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As  with  the  transformation  sequence,  the  steps  for  reduction  are  as  follows: 

1 .  Define  a  variable,  (P),  that  represents  the  number  of  DOF  of  the  reduced  system  model. 

The  only  caveat  is  that  the  selection  of  (P)  must  be  divisible  by  2.  The  resulting  state  space 
matrices  B  and  B  will  have  dimensions  (2P  x  2P).  This  requires  (P)  to  be 

^  Isys 

smaller  than  the  lowest  subsystem  DOF. 

2.  Reduce  the  size  of  modal  matrices,  T*  and  ,  to  dimensions  {2Ni  +  . . .  +  2Nm  x  2P)  by 

discarding  the  unnecessary  columns  of  the  modal  matrices.  Since  this  is  a  state-space 
formulation,  the  columns  associated  with  each  of  the  (2m)  vector  components  (i.e., 

{■i. {-i.  {'id*...  ■  “  "‘1“' 

reduction.  The  size  of  the  reduced  system  is  {2P  x  2P).  The  selection  matrices  from  T* 
and  T*  will  result  in  an  equal  number  of  modes  selected  from  the  original  eigen  matrices 
associated  with  each  vector  component.  The  reduced  transformation  matrices,  and 
,  will  thus  have  dimensions  {2N\  +  . . .  +  2Nm  x  2P). 

3.  Transform  the  system  matrices,  A  and  A  ,  using  the  reduced  transformation 

L  J  sys  L  ^  J  sys 

matrices,  [t*|  and  [t^|  •  This  results  in  the  reduced  diagonal  matrices  of  equation  45. 


r  tT 

T, 

.  J(2Px2y+---+2Ar.^) 

a; 

sys{2N-^-\ - h27V.^x2y  H - h2Af.^ ) 

t; 

'2y+---+2Ar.^x2P) 

b; 

sysi2Px2P 

»lT 

T, 

^J(2Px2y+---+2Ar.^) 

a; 

sys[2N^+---+2N^y.2N^+---+2N\ 

t; 

'2y +---+27V...X2P) 

K 

syshPx2P 

(45) 


4.  Solve  the  eigenvalue  problem  of  the  reduced  state-space  matrices,  and 

L  ^  Js^5(2Px2P) 

B  The  eigenvalues  will  result  in  {P)  complex  conjugate  pairs,  which  represent 

^J*as(2Px2Pj 

the  first  {P)  frequencies  of  the  system.  They  will  have  the  same  or  very  near  the  same 
value  as  the  corresponding  frequencies  of  the  full  system  model. 


3,2  Speed-Sweep  Vibration  Response 

The  determination  of  the  vibration  response  is  possible  using  the  principal  coordinates  as  well. 
Returning  to  equation  44,  one  can  assume  a  harmonic  response  to  the  forcing  function, 
{F(t)=^}.,.: 


such  that  the  equations  of  motion  become 
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(47) 


The  displacements  in  modal  coordinates  are 


Converting  back  to  physical  coordinates  yields  the  following  expression: 


(48) 


(49) 


where  is  the  vector  containing  the  amplitudes  of  the  response.  Substituting  the 

operational  speeds,  (ff),  into  equation  49  will  yield  speed  sweep  vibration  responses,  as  shown 
later  in  figures  3  and  4. 

Gyroscopic  effects  play  a  unique  role  in  speed-sweep  response.  The  modal  transformation 
matrices  are  determined  by  calculating  the  eigenvectors  of  the  system.  The  matrix,  [D],  contains 
the  system  gyroscopics,  which  is  a  function  of  rotational  speed.  Hence,  the  eigenvectors 
themselves  are  a  function  of  rotational  speed,  and  the  transformation  matrices  themselves  change 
with  rotational  speed.  Consequently,  for  a  full,  modal  speed-sweep  analysis,  the  transformation 
matrices  must  be  calculated  at  each  speed  interval  within  the  desired  range.  Additionally,  the 
frequency  response  can  be  calculated  by  fixing  the  operational  speed  in  the  [D]  matrix  and 
varying  the  frequency  {co)  in  equation  49. 

3.3  Time  Domain  Response 

If  desired,  it  is  also  possible  to  determine  the  time  domain  forced  response  using  equation  49. 

The  time  domain  response  acts  at  a  single  frequency.  The  shaft  operating  speed  is  simply 
substituted  into  equation  49  to  solve  for  the  amplitude  vector,  {Q*}eye.  The  time  domain 
response  is  simply  the  periodic  function: 

{q}  ={q1  e'"  (50) 

^  J  sys  ^  J  sys 
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4.  Application:  Geared  System  with  Damping 


4,1  Overview 

The  application  for  preliminary  validation  of  this  methodology  is  a  basic  model  of  a  spur  gear 
shaft  pair  similar  to  the  one  found  in  figure  1 .  Selected  because  of  the  ready  availability  of  the 
system  parameters,  it  appears  a  number  of  times  in  the  literature  (15,  16).  The  system  consists  of 
two  identical  spur  gears  resting  at  the  midpoint  of  identical  shafts,  each  connected  to  the  gearbox 
by  a  pair  of  rolling-element  bearings.  Since  the  two  substructures  are  the  same,  the  reduction  in 
the  complexity  of  equation  30  is  significant.  All  substructure  matrix  components  are  identical 
for  each  substructure.  Table  1  provides  the  parameters  of  the  two-shaft  system. 


Figure  1.  Spur  gear  test  stand  (75). 
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Table  1 .  Spur  gear  application  system  parameters. 


Material  Parameters 

Young’s  Modulus  (N/m^) 

2.03E11 

Shear  Modulus  {N/rn) 

8.0E10 

Poission’s  Ratio 

0.3 

Density  (kg/m^) 

7750 

Shaft  Parameters 

Length  (m) 

0.254 

Mass  (kg) 

1.96 

Rigidity  (N-in^) 

18576 

Outer  Diameter  (m) 

0.037 

Inner  Diameter  (m) 

0.01 

Gear  Parameters 

Mass  (kg) 

1.84 

Moment  of  Inertia  {kg-rr7) 

1.8E-3 

Polar  Moment  of  Inertia  (kg-m^) 

3.6E-3 

Diameter  (m) 

0.089 

Pressure  Angle  (deg) 

20 

Helieal  Angle  (deg) 

0.0 

Average  Mesh  Stiffness  (N/m) 

1.0E8 

Bearing  Parameters 

kxx,  kyy  (N-s/m) 

1.0E9 

Cxx  (N-s/m) 

500 

c,,  (N-s/m) 

700 

The  gear  mesh  is,  in  actuality,  a  periodic  (non-linear)  stiffness  due  to  the  alternating  number  of 
teeth  in  contact  throughout  the  mesh  cycle.  This  closed-form  solution  linearizes  the  equations  of 
motion  by  assuming  an  average  mesh  stiffness  value  (table  1).  Otherwise,  for  non-linear 
treatment,  numerical  integration  techniques  become  necessary  {4) 

The  full  finite  element  system  model  calculates  the  eigen  solution  of  the  system  and  the 
frequency  response  to  a  harmonic  forcing  function.  Each  shaft  consists  of  126  DOF,  or  252  for 
the  entire  system.  Since  the  state-space  formulation  doubles  the  size  of  the  model,  the  total 
number  of  DOF  is  504.  The  software  used  in  this  formulation  was  MATFAB®. 

Prior  to  any  modal  reduction,  the  validity  of  the  system  model  was  accessed  through  a 
comparison  of  natural  frequencies  at  low  speed  and  low  damping  to  those  previously  published 
in  the  literature  {16,  17).  As  depicted  in  figure  2,  the  model  provides  results  consistent  with 
previous  experiments.  The  x-axis  of  figure  2  contains  no  values.  It  is  a  set  of  bins  that 
categorizes  the  like  frequencies  of  the  three  sets  of  data. 
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Figure  2.  Comparison  with  previous  results  {16,17). 

4,2  Modal  Reduction  Comparisons 

Since  modal  analysis  is  an  approximation,  it  requires  further  validation  through  eomparison  with 
the  full  system  model  results.  The  first  validation  eompares  the  eigenvalues  of  the  system  model 
with  those  of  the  redueed  model.  Tables  2  and  3  provide  the  results  of  reduetions  to  (P)  =  30  and 
20  DOF,  respectively.  The  rotational  speed  (O)  for  this  analysis  is  1000  rpm  or  16.67  Hz.  For 
brevity,  only  the  first  10  frequeneies  are  provided. 


Table  2.  Modal  analysis  eigenvalue  comparison  (P=30). 


Modal  Synthesis 
(Hz) 

Full  FEA 
(Hz) 

Difference 

(%) 

1 

504.647 

504.635 

0.00% 

2 

684.261 

684.260 

0.00% 

3 

684.279 

684.279 

0.00% 

4 

684.296 

684.297 

0.00% 

5 

2091.562 

2090.630 

0.04% 

6 

2837.572 

2837.560 

0.00% 

7 

2837.572 

2837.560 

0.00% 

8 

2858.979 

2858.969 

0.00% 

9 

2858.979 

2858.969 

0.00% 

10 

6677.140 

6676.739 

0.01% 
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Table  3.  Modal  synthesis  eigenvalue  comparison  (P=20). 


Modal  Synthesis 

m) 

Full  FEA 
{Hz) 

Difference 

(%) 

1_ 

504.65204 

504.635 

0.00% 

1_ 

684.261 

684.260 

0.00% 

3 

684.279 

684.279 

0.00% 

4_ 

684.296 

684.297 

0.00% 

2091.950 

2090.630 

0.06% 

6_ 

2837.572 

2837.560 

0.00% 

1_ 

2837.572 

2837.560 

0.00% 

8_ 

2858.979 

2858.969 

0.00% 

2858.979 

2858.969 

0.00% 

10 

6677.140 

6676.739 

0.01% 

These  results  elearly  show  the  levels  of  aceuracy  attainable  by  modal  synthesis.  Furthermore, 
the  modal  simplifieations  also  produce  substantial  reductions  in  the  amount  of  computing  time 
required  compared  to  the  full  system  model;  62%  and  63%,  respectively. 

The  second  validation  method  occurs  by  comparing  the  reduced  system’s  speed-sweep  responses 
to  excitation.  The  responses  of  two  harmonic  forcing  functions  are  presented.  The  first  is  the 
gear-meshing  frequency.  The  gear-mesh  frequency  is  a  result  of  the  clatter  between  gears  as  the 
gear  teeth  come  in  and  out  of  contact  with  each  other  and  the  nonlinear  loading  of  the  gear  mesh. 
It  is  present  in  all  gears  and  for  spur  and  helical  gears,  it  is  defined  as  the  number  of  gear  teeth 
multiplied  by  the  shaft  rotational  seed.  The  key  impact  of  the  gear-mesh  frequency  is  that  it 
excites  higher  modes  of  vibration  at  lower  rotational  speeds. 

The  second  forcing  function  models  a  mass  imbalance  in  one  of  the  gears.  A  finite  imbalance 
exists  in  all  rotating  systems  and  is  a  function  of  the  rotational  speed  squared,  (ff)  .  Figures  3 
and  4  provide  the  speed-sweep  responses  for  the  gear-mesh  frequency  and  mass  imbalance 
excitations,  respectively,  using  a  speed  step,  (Aff),  of  50  Hz. 

The  response  curves  represent  the  lateral  displacement  of  each  shaft  at  its  midpoint.  Three 
curves  are  plotted  for  each  shaft,  providing  the  results  of  (1)  the  full  finite  element  model  -  “Full 
FEA,”  (2)  the  full  modal  reduction  -  “Modal  Synthesis,”  and  (3)  the  modal  reduction  using  a 
fixed  operating  speed  -  Modal  Fixed  Gyro.”  The  Modal  Synthesis  scenario  recalculates  the 
transformation  matrices  for  each  rotational  speed  interval.  The  Modal  Fixed  Gyro  scenario 
calculates  the  matrices  only  once  for  a  single  value  within  the  desired  speed  range  in  order  to 
reduce  the  number  of  burdensome  eigenvector  calculations.  It  then  uses  these  matrices  for  all 
responses  within  the  operating  range,  0-5000  Hz.  Figures  3  and  4  show  that  for  this  system,  the 
difference  between  these  two  approaches  is  negligible. 
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Shaft  1  Full  FEA  — B —  Shaft  1  Full  Modal  — —  Shaft  1  Fixed  Gyro - Shaft  2  Full  FEA  — ©—  Shaft  2  Full  Modal  — H —  Shaft  2  Fixed  Gyro 


Figure  3.  Gear-mesh  frequency  response  comparison  (0-5000  Fiz). 


Figure  4.  Imbalance  response  comparison  (0-5000  Hz). 
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If  necessary,  the  time  domain  response  can  also  be  determined.  This  might  be  desirable  if  the 
shape  of  one  of  the  shaft  orbits  is  important.  Figure  5  shows  four  views  of  the  combined 
horizontal  and  vertical  vibrations  for  a  single  shaft.  The  top-left  plot  simply  shows  the  circular 
motion  of  the  vibration  over  time.  The  top-right  chart  shows  the  shape  of  the  orbit.  The  bottom 
two  charts  simply  show  the  sinusoidal  motion  of  the  vibration  in  the  horizontal  (left)  and  vertical 
(right)  planes  of  motion. 


Figure  5.  Time  response  of  damped  gear  shaft  due  to  excitation  near  resonance. 


5.  Discussion  -  Modal  Synthesis  Practicality  in  Geared  System  Application 


For  all  its  accuracy,  this  experiment  does  raise  some  questions  regarding  the  utility  of  such 
procedure  in  terms  of  practicality  and  efficiency.  These  questions  fall  into  two  categories; 
limitations  on  modal  reduction  and  the  impact  of  the  gyroscopic  effect. 
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5,1  Limitations  on  Modal  Reduction 


Once  again,  the  modal  synthesis  procedure  can  produce  a  signifieant  modal  reduction  without 
loss  of  accuracy.  However,  in  the  speed-sweep  response,  the  gear-mesh  forcing  function  will 
limit  the  amount  of  reduction  available.  At  lower  speeds,  the  gear-mesh  frequency  will  excite 
higher  modes  of  vibration.  Consequently,  for  the  modal  synthesis  transformation  to  capture 
these  higher  modes,  they  must  be  included  in  the  modal  reduetion.  An  example  is  presented  in 
figure  6. 


Rotational  Speed  (H)  (Hz) 


Shaft  1  Full  FEA  - 


-  Shaft  1  Full  Modal  - 


-  Shaft  1  Fixed  Gyro  - 


-Shaft  2  Full  FEA  - 


-  Shaft  2  Full  Modal  - 


-  Shaft  2  Fixed  Gyro 


Figure  6.  Gear-mesh  frequeney  limitations  on  modal  reduetion. 

The  total  number  of  DOF  in  the  geared,  two-shaft  model  is  252.  This  increases  to  504  DOF  in 
the  state-spaee.  Figure  6  illustrates  a  modal  reduction  to  30  DOF  (60  DOF  in  the  state-space)  in 
terms  of  the  gear-mesh  forcing  function  response.  The  reduced  model  has  failed  to  accurately 
portray  all  resonance  peaks  within  the  operating  range,  0-5000  Hz.  To  capture  the  remaining 
peaks,  the  modal  DOF  must  increase  to  aecount  for  excitation  of  the  higher  modes.  For  the 
application  system,  the  required  number  of  degrees  of  freedom  to  produce  an  accurate  response 
is  168,  or  in  the  state-space,  336.  Thus,  an  intended  88%  DOF  reduetion  will  result  in  a 
reduction  of  only  33%  if  the  gear-meshing  frequeney  excitation  is  important. 

The  implication  for  geared  system  analysis  is  that  there  must  be  some  a  priori  knowledge  as  to 
the  number  of  higher  modes  excited  within  the  desired  operating  speed  range.  An  increase  in  the 
number  of  gear  teeth  or  the  operating  range  itself  will  increase  the  number  of  required  modes. 
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The  amount  of  modal  reduction  available  is  therefore  inversely  proportional  to  the  number  of 
gear  teeth  and  shaft  operating  speed. 

5.2  Gyroscopic  Impacts 

The  presence  of  gyroscopics  complicates  the  analysis  of  any  system,  and  modal  synthesis  is  no 
different.  As  previously  discussed,  a  full  modal  speed-sweep  response  analysis  requires  the 
transformation  matrices  to  be  calculated  at  each  speed  interval,  which  greatly  adds  to  computer 
processing  time.  Nevertheless,  figures  3  and  4  showed  that  one  could  reach  a  good 
approximation  by  “fixing”  the  gyroscopics  and  calculating  the  transformation  matrices  at  one 
speed  only.  However,  the  ability  to  fix  the  gyroscopics  depends  upon  the  system’s  sensitivity  to 
gyroscopic  effects. 

5.3  Limitations  to  Practicality 

These  two  factors  may  adversely  affect  the  utility  of  the  methodology  by  increasing  the  computer 
processing  time  of  a  procedure  originally  designed  to  reduce  it.  A  comparison  of  computing 
time  is  presented  in  table  4,  illustrating  the  gyroscopic  effects  on  processing  time.  The  full  252- 
DOF  system  model  is  the  comparison  standard  with  a  computer  processing  ratio  of  one.  In  the 
Full  Modal  Synthesis  models,  where  the  transformation  matrices  are  recalculated  for  each 
operating  speed  within  the  spectrum,  the  168  and  30  DOF  models  take  8.7  and  5.9  times  longer 
to  process,  respectively,  than  the  full  finite  element  model.  The  Modal  Fixed  Gyro  scenarios 
reduce  those  times  substantially.  In  the  case  of  the  30  DOF  model,  the  processing  time  is  20%  of 
the  full  system  model.  However,  as  pointed  out,  these  two  reduction  models  do  not  produce  an 
accurate  gear-mesh  frequency  response. 


Table  4.  Vibration  response  processing  time  ratio  comparison 


Degrees  of  Freedom 

252 

168 

30 

Full  Modal  Synthesis  Scenario 

1 

8.7 

5.9 

Modal  Fixed  Gyro  Scenario 

1 

1.3 

0.2 

Consequently,  while  this  experiment  validates  the  modal  synthesis  procedure  and  results  in 
computer  processing  time  savings  during  eigen  analysis,  the  cutoff  frequency  or  limit  of  modal 
reduction  is  determined  by  the  type  of  excitation  analyzed.  Therefore,  it  may  or  may  not  be 
particularly  practical  for  use  in  vibration  response  calculations. 
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6.  Conclusions 


This  report  has  presented  a  method  for  performing  modal  synthesis  on  a  geared  rotor  dynamie 
system  undergoing  general  damping  and  gyroseopic  effects.  The  modal  analysis  produces  the 
eigen  solutions,  frequency  response,  and  time  domain  response  of  the  system.  Validation  of  the 
reduced  model  occurred  by  comparing  the  eigen-solutions  and  frequency  responses  of  the 
reduced  system  with  the  original  model.  The  results,  in  almost  all  instances,  matched  exactly. 

In  the  speed-sweep  response,  the  modal  synthesis  procedure  resulted  in  limitations  due  to  the 
limit  of  modal  reduction  imposed  by  gear-mesh  frequency  considerations  and  the  amount  of 
computer  processing  time  required.  The  gear-mesh  frequency  excited  higher  modes  of  vibration, 
which  required  a  higher  DOF  model  to  maintain  accuracy.  The  introduction  of  gyroscopic 
effects  increased  processing  time,  substantially  in  some  cases.  These  two  limitations  call  into 
question  the  applicability  and  practicality  of  using  the  modal  synthesis  method  for  geared, 
rotating  systems. 

Finally,  an  important  point  is  that  this  was  only  one  system  analysis;  other  systems  may  behave 
differently.  For  example,  fixing  the  gyroscopic  speed  value  in  the  transformation  matrix  of  this 
scenario  was  acceptable.  Other  systems  may  be  more  susceptible  to  gyroscopic  effects,  where 
fixing  the  speed  would  be  unacceptable.  It  also  may  be  that  in  other  systems,  the  ratio  of 
computer  processing  time  is  substantially  lower  for  response  calculations,  making  it  an  attractive 
alternative  to  the  full  finite  element  model,  especially  when  the  full  model  contains  thousands  or 
millions  of  DOF. 


7.  Recommendations 


Recommended  further  research  falls  into  three  categories.  The  first  deals  with  application  to 
other  systems.  This  procedure  worked  very  well  for  a  simple,  symmetric,  two-shaft  system. 
More  complex  systems  might  behave  differently,  especially  when  gyroscopic  sensitivity  is 
considered.  Second,  the  extension  of  this  method  in  complex  notation  to  capture  forward  and 
backward  whirl  might  also  be  of  interest.  Finally,  this  method  should  be  further  considered  for 
its  potential  to  other  applications,  including  system  identification  and  design  optimization 
methods. 
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List  of  Symbols 


j 

K 


m. 


'I 


N 

P 


{^} 

{4 

{-I 


general  damping  matrix 

bearing  damping  values  inx  and  j-directions,  respectively 

combined  damping/gyroscopic  matrix 

diagonal  damping  matrix  component 

excitation  force  vector 

gyroscopic  matrix 

identity  matrix 

imaginary  component,  P~i 
stiffness  matrix 

diagonal  stiffness  matrix  component 

bearing  stiffness  values  inx  andy-directions,  respectively 

transformation  sub-matrix  designation 

consistent  mass  matrix 

diagonal  mass  matrix  component 

number  of  degrees  of  freedom 

number  of  reduced  degrees  of  freedom 

amplitude  of  physical  coordinate  vector  of  general  equation 

of  motion 

physical  coordinate  vector 

first  derivative  of  physical  coordinate  vector 

first  derivative  of  principal  coordinate  vector 
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a 


P 


c 


A 


LU 


principal  coordinate  vector 
original  state-space  matrix 
transformed  state-spaee  matrix 
transformation  matrix 

shaft  rotational  speed 
real  eigenvalue  eomponent 
imaginary  eigenvalue  eomponent 

diagonalizing  matrix 

damping  ratio 

diagonal  eigenvalue  matrix 

eomplex  eigenvalue 
elimination  matrix 

redueed  system  modal  matrix 
original  modal  matrix 
eigenvector 
system  frequeney 


Subseripts 

D 

diagonal 

Dm 

transformed  mesh 

1 

left 

m 

number  of  substruetures  of  global  system 

mesh 

gear-mesh 

n 

natural  frequeney 

r 

right 

sys 

global  system  to  inelude  eoupling  effects 
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Superscripts 

T  transpose 

'  first-order  vector,  scaled  or  normalized 

*  global  system  excluding  coupling  effects 
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